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We obtain a good analytic fit to the joint Bjorken-x and Q 2 dependences of ZEUS data on the 
deep inelastic structure function F2(x,Q 2 ). At fixed virtuality Q 2 , as we showed previously, our 
expression is an expansion in powers of ln(l/a;) that satisfies the Froissart bound. Here we show 
that for each x, the Q 2 dependence of the data is well described by an expansion in powers of InQ 2 . 
The resulting analytic expression allows us to predict the logarithmic derivatives (c?"i 7 ' : f / (d In Q 2 ) n ) x 
for n = 1,2 and to compare the results successfully with other data. We extrapolate the proton 
structure function F%(x, Q 2 ) to the very large Q 2 and the very small x regions that are inaccessible 
to present day experiments and contrast our expectations with those of conventional global fits of 
parton distribution functions. 
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Introduction. The ability to predict cross sections at 
very high energies, whether at the CERN Large Hadron 
Collider or in ultra high energy cosmic ray interactions, 
depends critically on the reliability of extrapolations from 
current measurements into regions of much greater vir- 
tuality (Q 2 ) of the elementary scattering processes, and 
to much smaller values of the fractional longitudinal 
momentum x carried by the parton constituents of the 
hadrons. Most high energy predictions are expressed in 
terms of convolutions of elementary hard-scattering cross 
sections with parton distribution functions (PDFs). The 
quantitative reliability of these predictions relies on the 
x and Q 2 dependences embodied in the universal par- 
ton distribution functions extracted from global analy- 
ses in perturbative quantum chromodynamics (QCD) of 
data at lower energies. Important in these global analy- 
ses are data from deep-inelastic lepton scattering (DIS) 
and other reactions, but equally critical are the analytic 
functional forms assumed for the x dependence of parton 
distribution functions at small x. 

In an earlier paper [l[ , we analyze the x dependence of 
the DIS proton structure functions -Ff (x, Q 2 ). We begin 
with the assumption that the x dependence at extremely 
small x should manifest a behavior consistent with sat- 
uration of the Froissart bound on hadronic total cross 
sections [2j, as is satisfied by data on yp, n p, and pp 
and pp interactions |3| . This bound [2j , derived from an- 
alyticity and unitarity, demands that Ff (x, Q 2 ) grow no 
more rapidly than ln 2 (l/x) at very small x for all val- 
ues of Q 2 . Over the ranges of x and Q 2 for which DIS 
data are available, we then show that a very good fit 
to the x dependence of ZEUS data ;4j is obtained for 
x < xp = 0.09 and — 3> m 2 with the expression 
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Our fits to DIS data [4| at 24 values of Q 2 cover the wide 
range 0.11 < Q 2 < 1200 GeV 2 . The value x P = 0.09 is 
a scaling point l] such that the curves for all Q 2 pass 
through the point x = xp, at which F^xp, Q 2 ) = Fp ~ 
0.41, constraining all of the fits. 

In this note, we extend our analysis by making a joint 
fit of both the x and Q 2 dependences of the ZEUS [3| 
data on F® (x, Q 2 ). The analytic expression we derive for 
the x and Q 2 dependences allows us to compute the loga- 
rithmic partial derivatives [d n F^{x, Q 2 ) / (dh\Q 2 ) n ) for 
n = 1,2. We obtain excellent agreement when comparing 
our predictions for the first derivative (n = 1) with HI 
data [H}. We offer predictions for the second derivative 
(n = 2). Only 8 parameters — two of which are the scal- 
ing value Fp = 0.41 at the scaling point xp — 0.9 — are 
needed to fit the joint x and Q 2 dependences. Our ex- 
pression allows us to extrapolate to very low values of x, 
well beyond the experimental range presently accessible. 
We obtain cross sections for ultra high energy cosmic ray 
neutrino reactions that are a factor of ~ 5 smaller than 
those based on extrapolations of conventional parton dis- 
tribution functions. 

Joint Fit. In Eq. (jTJ) , Ff (x, Q 2 ) is written as a sum of 
terms that are factorizable as functions of Q 2 times func- 
tions of x that satisfy the Froissart bound. In this paper 
we discuss our fit for the functions A{Q 2 ) and B(Q 2 ). 

To parameterize the dependence on Q 2 at fixed x, we 
assume an expansion in powers of InQ 2 , generally con- 
sistent with and motivated by the InQ 2 variation ex- 
pected in QCD. We note, moreover, that the HI collabo- 
ration H determined that, for fixed x, the Q 2 dependence 
of Fg(x,Q 2 ) is reproduced by the form Fg(x,Q 2 ) — 
ao(x) + ai(x) ln(<3 2 ) + 02 (x) ln 2 (Q 2 ). We therefore ex- 
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pand the functions A(Q 2 ) and B{Q 2 ) as 
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terminating these phenomenological expansions at the 
quadratic level. 

We fit simultaneously the x dependence and the Q 2 de- 
pendence of the data (Q 2 is expressed in GeV 2 through- 
out). We determine the 6 real constants a_, Oj., 0,2, bo, b\ 
and b 2 in Eq. ([2]) using the Sieve algorithm [||, by mini- 
mizing the squared Lorentzian, 
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where x 2 (a;x) = V\ =1 A%f (ar»; a), Ax 2 (^;a) = 
— yi(xj)] /(Tj) , a is the parameter space vec- 
tor, and a) is the theoretical value of the measured 
2/j at x_, with measurement error <x,-. Using a A% 2 max cut 
of 6, we find (see Table[T| a final corrected x 2 /d.f.=1.09, 
for 169 degrees of freedom (d.f.), a reasonable fit. The 
data used are 24 ZEUS data sets, with Q 2 =0.11, 0.25, 
0.65, 2.7, 3.5, 4.5, 6.5, 8.5, 10, 12, 15, 18, 22, 27, 35, 
45, 70, 90, 120, 200, 250, 450, 800 and 1200 GeV 2 . The 
Sieve algorithm eliminated 8 outlier points (from a total 
of 183) which had a \ 2 contribution of 63.5. 

The quality of our fit to the x and Q 2 dependences of 
the data for x < xp, — 3> to 2 is shown in Fig. [T] a 
representative plot of 13 of the data sets. It shows that 
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FIG. 1: Fits to the proton structure function data, F^ix, Q ) 
vs. x, for 13 values of Q 2 . The data are from the ZEUS 
collaboration Q. The curves show 13 of our 28 global 
fits whose parameters are given in Table __. The vertical 
and horizontal straight lines intersect at the scaling point 
x P = 0.09, Ff (z P ) = 0.41. 

the fit is excellent over the large Q 2 and x ranges of the 
ZEUS Q data, fitting data equally well at x ~ 10~ 6 for 



Q 2 = 0.11 GeV 2 as at x - 10~ 2 for Q 2 = 1200 GeV 2 , 
even with constraint that all curves must pass through 
the common scaling point xp — 0.09, Fp — 0.41. The 
6 fit parameters are given in Table [Q along with their 
statistical errors. 

TABLE I: Results of a 6-parameter fit to F% (x, Q 2 ) structure 
function dataQ using the x and Q 2 behaviors of Eq. (_[]) and 
Eq. ([__), with Q 2 in GeV 2 . The renormalized Xmin P er degree 
of freedom, taking into account the effects of the Ax 2 max = 6 
cut [_]], is given in the row labeled 1Z x Xmin/dT The errors in 
the fitted parameters are multiplied by the appropriate r X 2[l]. 



Parameters 


Values 


a 


-5.381 x 10 -2 ±2.17 x 10 -3 


0,1 


2.034 x 10~ 2 ± 1.19 x 10~ 3 


ai 


4.999 x 10~ 4 ±2.23 x 10" 4 


bo 


9.955 x 10" 3 ±3.09 x 10~ 4 


61 


3.810 x 10" 3 ± 1.73 x 10" 4 


&2 


9.923 x 10~ 4 ±2.85 x 10" 5 


Xmin 


165.99 


^ X Xmin 


184.2 


d.f. 


169 


11 X Xmm/d-f. 


1.09 



Evaluation of dF% (x,Q 2 )/d\n(Q 2 ). Differentiating 
Eq. Ul with respect to ln(Q 2 ), we obtain 

dF 2x( x ,Q 2 ) 
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an expression valid for x < xp and Q 2 /x ^> m 2 . We 
show a plot of dFl x (x,Q 2 )/d\n(Q 2 ) in Fig. [2 for a set 
of values of x, compared to the values measured by the 
HI collaboration [5j. We emphasize that the theoretical 
values have been constrained by ZEUS [3] data alone, 
and that they are a prediction of the HI results, not a fit 
to these data. We see from Fig. __] that our predictions 
based on the ZEUS data are in fine agreement with the 
normalization and slope of the HI results. 

Curvature. The curvature, defined as one-half of the 
second derivative, is given by 
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The form of Eq. |5]) indicates that our curvature is inde- 
pendent of Q 2 , a consequence of the fact that we truncate 
the expansions in Eq. @ at the quadratic level. With our 
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FIG. 2: A plot of the derivative dF$ x (x, Q 2 )/d ln(Q 2 ) vs. Q 2 , 
in GeV 2 , for selected values of x, compared to data from the 
HI collaboration The exterior lines for x = 0.00008 are 
the error bands associated with the parameter uncertainties 
of the coefficients of Table [T] 



parameterization, the curvature grows like ln 2 (l/a;) as x 
decreases. The signs and magnitudes of a 2 and hi de- 
termine the sign (positive/negative) of the curvature. In 
our case, the curvature is positive, and it increases as x 
decreases, features that are also true in next-to-leading 
order perturbative QCD [7|. We note, however, that we 
do not impose or employ QCD evolution in our work. 
Our predictions are based entirely on our fit to data on 
F%(x,Q 2 ). The results of our calculation of curvature 
are in good agreement with the data shown in Ref. Q. 
We remark that curvature is defined somewhat differ- 
ently in Ref. Q as the second derivative with respect to 
log 10 (l + Q 2 /Qo); instead of with respect to ln(Q 2 ). 

Extrapolation to Very Small x. In Fig. [31 we present 
our calculation of F^x^Q 2 ) as a function of x for the 
choices Q 2 = 25 GeV 2 and Q 2 = M w , where M w 
is the mass of the intermediate W boson. The scale 
Q = M\y is of interest at hadron colliders where it char- 
acterizes electroweak processes. It is also the relevant 
scale in charged-current high energy neutrino interac- 
tions [9j where the W boson propagator limits momen- 



tum transfers to Q 2 ~ M 2 ^. 

We contrast our expectations with evaluations of 
F2(x,Q 2 ) based on the CTEQ6.5 set of parton distri- 
bution functions (loj . In our case, the uncertainty bands 
represent a ±3 standard deviation variation of our pa- 
rameters, whereas in the CTEQ6.5 case the bands are 
obtained from the 40 eigenvector sets that encapsulate 
the uncertainties of their PDFs. 

Inclusive W boson production at hadron colliders 
serves as an independent probe of the x dependence of 
quark and antiquark densities at Q — J%, sensitive to 
values of x ~ (Mw / \/s)e~ v , where s is the square of the 
center-of-mass energy and y is the rapidity of the W bo- 
son [ll|. At the Fermilab Tevatron, with y/s = 1.96 TeV, 
values of x down to 2 x 10 ~ 3 are probed for y ~ 3, well 
within the range in which our expectation and that of 
CTEQ6.5 show agreement in Fig. [3J 

The magnitude and x dependence of the CTEQ6.5 
and our calculations agree quite well over the range 
10~ 3 < x < 0.1. Both approaches also show the same 
dependence on Q 2 in this region of x. The agreement 
is expected since both fit data that are limited to this 
range of x at large Q 2 . The agreement also shows that 
the logarithmic expansion we use to describe x depen- 
dence and the inverse power behavior of the CTEQ form 
cannot be distinguished numerically over the finite range 
10~ 3 < x < 0.1. However, the two expectations clearly 
diverge considerably when extrapolated to values of x as 
low as 10~ 8 . In the absence of new physics effects in 
the small- a; region, the saturation of the Froissart bound 
embodied in our fit to data at available energies yields a 
robust extrapolation into the region of very small x. 

In the parton model, the decomposition of the struc- 
ture function F%{x, Q 2 ) at very small x is dominated by 
the sea quark q(x,Q 2 ) and sea antiquark q(x,Q 2 ) den- 
sities. Although we do not decompose our F^ into par- 
ton distributions, this dominance allows us to conclude 
from Fig. [3] that our sea quark (and antiquark) distribu- 
tions will be ~ 5 times smaller than those in CTEQ6.5 at 
x ~ 10~ 8 and Q 2 ~ M w . In ultra high energy charged- 
current neutrino interactions x ~ M w /2mE u , and 
the cross section on (isoscalar) nucleons is proportional 
to the sum of the up and down antiquark distributions, 
u(x,M w ) + d(x,M w ). From Fig. [31 we see that the 
expected cross section at x ~ 10~ 8 will be ~ 5 times 
smaller than that based on CTEQ6.5 [lil] . serving to de- 
fine a range of uncertainty for estimated rates, with our 
expectation coming in on the low side. 

Summary. The Bjorken-x dependence of the DIS pro- 
ton structure function F^x, Q 2 ) measured by the ZEUS 
collaboration is consistent with a ln 2 (l/x) dependence 
at small values of x, compatible with saturation of the 
Froissart bound at each value of Q 2 . We parameterize 
successfully the joint x and Q 2 dependences of F2 for 
x < xp ~ 0.09 and — 3> m 2 , using the compact fac- 
torized expression in Eq. JT]), with the Q 2 variation ex- 
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FIG. 3: A plot of F?(x,Q 2 ) vs. x at Q 2 = 25 GeV 2 and 
Q 2 = M|^, where Mw is the mass of the W boson, along with 
results based on the CTEQ6.5M parton distribution func- 
tions fi~ot | . 

pressed in Eq. @ . Our analytic expression has only 6 pa- 
rameters (plus the scaling point xp and the value F^xp) 
at the scaling point). We compute the first and second 
derivatives of F^{x, Q 2 ) with respect to InQ 2 at small x. 
Our predictions of these quantities are in good agreement 
with the measurements by the HI collaboration We 
extrapolate our expression for F%{x, Q 2 ) down to the very 
small value x — 10~ 8 and compare our expectations to 
those based on the CTEQ6.5M parton distributions [lbT |. 

Under the assumption that the Froissart bound applies 
to the virtual photon cross section cr(7*p), a ln 2 (l/a;) be- 
havior is as singular as is allowed for the very small x be- 
havior of F% (x, Q 2 )- However, it is difficult to reconcile a 
ln 2 (l/x) behavior at very small x for F%{x, Q 2 ) [and for 
the gluon distribution g(x, Q 2 )] at all Q 2 with the Dok- 
shitzer Gribov Lipatov Altarelli Parisi (DGLAP) evolu- 
tion equation [l3| at next-to-leading order in QCD, ow- 
ing to the singular nature of parton splitting functions at 
small x. Global analyses of parton distribution functions 
based on DGLAP evolution, such as CTEQ6.5, begin 
with the assumption of an inverse power behavior for the 
small x dependences of the quark, antiquark, and gluon 
densities, behavior that is more singular than is allowed 
by the Froissart bound. The assumed inverse-power be- 
havior leads to the very different expectations shown in 
Fig. [3j where they are seen to diverge for x <, 10~ 3 . To 
the extent that the ln 2 (l/a;) behavior is preferable the- 
oretically, we question the reliability at very small x of 
parton distribution functions based on an assumed in- 
verse power behavior. 

The ln 2 (l/x) behavior that we show is consistent with 
the DIS data may be the signal for the onset of a new 
regime, the physics of saturation or gluon recombination 
processes g + g — > g at high parton densities. Thus, 



new theoretical effort is warranted to devise a QCD evo- 
lution framework compatible with ln 2 (l/x) behavior of 
parton densities at very small x, and experimental pro- 
grams should be pursued to measure the x and Q 2 vari- 
ations of structure functions at much smaller values of x 
than are currently explored. 

Our next goals include a reanalysis of all available data 
for F%{x,Q 2 ) and dF$ x (x, Q 2 )/dln{Q 2 ) for x < x P , in 
ep, \xp, and vp deep inelastic scattering, in order to ob- 
tain a joint fit to both the x and Q 2 dependences, con- 
strained by the Froissart bound and the scaling point. 
This work should allow us to make more accurate pre- 
dictions of the proton structure function at very small x 
and very large Q 2 , regions beyond the reach of existing 
accelerators. We also will investigate the domain of com- 
patibility in Bjorken-x of a In 2 \ jx behavior of F%(x, Q 2 ) 
at very small x with quark and gluon distribution func- 
tions obtained in the standard fashion, e.g., in Ref. [loj ]. 
with DGLAP evolution and an assumed inverse power 
behavior of PDFs. 
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